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Binding energy of shallow donors in a quantum well in the presence of a tilted 

magnetic field. 
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We present results of variational calculations of the binding energy of a neutral donor in a quantum 
well (QW) in the presence of a magnetic field tilted relative to the QW plane. Assuming that the 
donor is located in the center of the QW, we perform calculations for parameters typical of a 
II- VI wide-gap semiconductor heterostructure, using as an example the case of a rectangular CdTe 
quantum well with CdMgTe barriers. We present the dependence of the binding energy of a neutral 
donor on the tilt angle and on the magnitude of the applied magnetic filed. As a key result, we 
j /' show that measurement of the binding energy of a donor at two angles of the magnetic field with 

respect to the quantum well plane can be used to unambiguously determined the conduction band 
offset of the materials building up heterostructure. 
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I. INTRODUCTION 

Much work has already been done on calculating the energies and wavefunctions of electronic states in semiconductor 
quantum wells (QWs) in the presence of an applied magnetic field. Most of these theoretical studies have focused on 
the binding energy of the neutral donor (D°) 0, Q> IE Gl > charged donor 0,0, El; neutral exciton 0, IE El El > an d 

S charged exciton (trion) H EE EJ EE 
states as a function of magnetic field. The vast majority of these calculations 
treat the most straightforward case, where the magnetic field is applied parallel to the direction of growth of the QW. 
j Experimental and theoretical results show that in this geometry the binding energy of the above complexes increases 
(—] . with increasing magnetic field |9j, Il3| . In spite of considerable progress in this area, little attention has been paid to 
the dependence of the binding energy on the tilt angle of the field relative to the QW plane 0, fljl Il6| . 

The objective of the present work is to determine the dependence of the binding energy of a neutral donor on the tilt 
angle 9 between two limiting geometries, the first geometry (denoted below as Case 7) corresponding to the magnetic 
t-H , field B aligned along the growth axis of the QW (designated as the z-direction); and the second limit (denoted as 
Case 77) corresponding to external magnetic field applied in the plane of the QW. In our notation described below, 
\ in Case 7 we define the tilt angle as 9 = 0°, and in Case 77 as 9 = 90°, see Fig. ^ for details. 

It is well established that the binding energy of different electronic complexes stemming from the Coulomb inter- 
action increases as the dimensionality of the quantum structure decreases, i.e., as we progress from quasi-two- to 
i quasi-one- and eventually to quasi-zero-dimensional quantum structures . An external magnetic field localizes the 
charged particles in the plane perpendicular to B in the form of its cyclotron motion, while the particle can move 
freely in the direction of the applied field, constituting in effect one-dimensional localization EH . One should note 
that in this case the density of states also has the character of a one dimensional system, manifesting itself as peaks 
at the Landau level positions. For an electron subjected simultaneously to the potential of the QW and of an external 
magnetic field, 'total' localization of a particle is different in the two limiting cases defined above. In Case 7, the 
combined action of QW confinement and of magnetic localization have different directions, which then manifests itself 
as quasi-zero-dimensional localization. In Case 71, the QW and the magnetic confinements have the same direction, so 
that the electron retains its quasi-one-dimensional character associated with the magnetic confinement. This implies 
O ■ that the binding energy of a donor should be larger in Case 7 than in Case 77. When the tilt angle 9 increases, we 
can then say that the dimensionality of an electron is between quasi-zero and quasi-one, and we expect the binding 
energy of D° to be a monotonia function of the tilt angle 9. We expect analogous (monotonic) behavior as a function 
of the tilt angle for the binding energies of neutral excitons (composed of two charged particles) or trions (composed 
of three charged particles). These predictions are in contradiction to the results obtained for the binding energy of 
neutral donors in a GaAs quantum well presented in Ref. 

The present paper is organized as follows. We start with a model Hamiltonian and a class of trial wave functions. 
We then present results of variational calculations for the CdTe QW with CdMgTe barriers. These include calculations 
of the binding energy Ef, as a function of the tilt angle 9, as a function of the magnetic field B, and as a function 
of the quantum well barrier height V e . We show, finally, that measurement of the binding energy of D° in the two 
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FIG. 1: Proposed experimental geometry. CdTe QW is grown in the z-direction: L is the QW width and V e is the QW height. 
The donors are incorporated into the system by 5-doping technique only in the center of the well. The edge of the conduction 
band is sketched on the right hand side. The magnetic field B lies in the XZ plane. Case I corresponds to B — Be z and Case 
77 corresponds to B — Be x . 



limiting geometries (Case I and Case II on Fig. \Q can be used as a tool for determining the conduction band offset 
at the well/barrier interface. 



II. MODEL 



The Hamiltonian of a shallow donor embedded in a symmetric square quantum well is modeled by the following 
Hamiltonian 

H = + V^(z) -/--^r- » B 9lBs. (1) 

2m* 47ree |f _ R Q \ 

The first part of the Hamiltonian is the kinetic energy of a delocalized conduction electron (where e is the electron 
charge, and m* is its effective mass) in presence of a tilted magnetic field B = B(sm(6), 0, cos(0)) lying in the XZ 
plane, see Figf^ We have chosen an asymmetric gauge for vector potential A(f) = B(0,xcos(8) — zsin(9), 0). For 
8 = 0° (Case 7), the magnetic field B is parallel to the OZ axis; and for 8 = 90° (Case II) it lies in the XY plane (i. 
e., in the plane of the QW), see Figf^ The profile of the potential energy of the QW is described by the second term 
in Eq. JTJ): 

where L is the width of the QW centered at z = 0, and V e is the height of its barrier. The energy scale is chosen 
by defining the conduction band edge of the barriers as zero. The third term in Eq. JIJ is the Coulomb energy of a 
shallow donor located at point Rq. We assumed that the donor center is located at the center of the QW, so we can 
set Ro = without losing physical generality. The last expression in Eq. (JIJ is the Zeeman Hamiltonian, in which 
g* is the effective g- factor of conduction electrons. We find it inconvenient to rotate the Hamiltonian, Eq. (JTJ, to a 
new coordinate system in which the kinetic energy has a simpler form than in the original system, because then the 
functional form of the potential profile, Eq. J2J), couples two variables (x and z in our notation) in a non-trivial way 

0E3. 

So far the problem of the Hamiltonian of a donor in a QW has not been solved analytically (even the case of a free 
electron in a QW in a tilted magnetic field remains analytically unsolved [ 19i \'2(l |21| ) . so that in our work we have 
used a variational approach. We propose the following form of the trial wave function of a two-component spinor, 

Ni N 2 N 3 

*± = f±(r)-X± ^EEE^^^i^H^' ( 3 ) 

i=0 j=0 k=0 
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where 4n are a one-dimensional harmonic oscillator functions (Gaussian functions) and x± are spin states. Note 
that, Gaussian trial wave functions have been successfully used in Ref. [l|, |l0| for calculating donor and trion states, 
respectively. The nonlinear variational parameter a (the scaling parameter), and the linear variational parameters 
^ijk were determined using the Ritz variational method. In Eq. the number of the basis functions has to be finite, 
and in this connection we have checked that Ai = A2 = A3 = 10 are sufficient to ensure that the results do not 
depend on the cut-off of the number of basis functions. 

The orbital part of the total wave function of a donor Eq. J2J), is denoted by f±(r), X± being the spin part. It 
is easy to show by direct substitution that the two spinors, \* — ( cos (^/2), sin(9/2)) and x = (—sin(6/2), cos(9/2)), 
solve the Schrodinger equation that contains the Hamiltonian given by Eq. Then the orbital part f±(r) of the 
spinor function i&± satisfies the following eigen-equation: 



Additionally, as seen from the above equation, f±(r) has the same functional form for both spin configurations x±- 



III. RESULTS AND DISCUSSION 



In order to demonstrate the tilt angle dependence of the binding energy of D , we performed calculations for 
two barriers heights, V e =200 meV and 20 meV. These two choices of barriers height correspond to 19% and 2% 
content of Manganese in the barriers, respectively. Additionally, for each V e we chose L=100 A and L=300 A; and 
B=0, 1, 4, 9, and 16 T, corresponding to magnetic lengths A c =oo, 256, 130, 85, and 65 A. For specificity we used 
CdTe materials parameters: m*=0.1 of free electron mass for both the QW and the barrier, and dielectric constant 
e=10.4, which give the characteristic Coulomb scales: i?y* = 12.6 meV and a|j=55 A. 

The donor binding energy Eb is obtained as a difference between the ground state energies of the free electron and of 
the donor (E±) with the same electron spin configuration [18j. This definition implies that the Zeeman Hamiltonian 
does not contribute to the binding energy of D°. We will be interested mainly in the variation of Eb with the 
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FIG. 2: Donor binding energy Eb in a rectangular CdTe/Cdo.8iMgo.i9Te quantum well as a function of tilt angle 8 of external 
magnetic field B for different magnetic field values (different symbols), calculated for two QW widths L. Full lines correspond 
to L—100 A, and dashed lines to L=300 A. Lines without symbols correspond to B=0 T. 

inclination 9 and the magnitude of the magnetic field, Eb — Eb(9, B). It must be noted that Eb also depends on other 
parameters, e. g., the barrier height, but these are kept constant in each variational process. 

In Fig. [21 we present the binding energy of the donor ground state as a function of tilt angle 9 for V e =200 meV at 
different magnetic fields, as well as for two different QW widths L=100 and 300 A. First we discuss the L=100 A case, 
represented by full lines in Fig. [21 
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At 73=0 T, the binding energy of the neutral donor is about 1.7 7?y*=21.4 meV, and is increased by 70% compared 
to the binding energy of a donor in three dimensions (in this case L ~ as). Next, for 73 ^ T, the binding 
energy is a monotonically decreasing function of the tilt angle: at a given magnetic field B = const, the binding 
energy is highest for 9 = 0°, and decreases for 9 > 0°. Suspection of Fig. [21 shows that for 9 = 0° the difference 
E b (B = \&T)-E b (B = 0T) is 0.6 Ry*=8 meV, while for 9 = 90° it is only 0.05 Ry*=0.6 meV, see also Fig. ©. These 
totaly different values in the two limiting field orientations (case 7 and case 77) are related to the fact that the QW 
width is smaller then (or comparable to) the characteristic magnetic length A c at fields up to 16 T. If the magnetic field 
is applied along the z-direction (6 = 0°), the electron is localized in x- and y-directions by the external magnetic field, 
as discussed at the outset (the bigger the field, the larger the magnetic localization). When this effect is combined 
with QW confinement, the electron becomes localized in all three direction. As the magnetic field is changing from 
T to 16 T, the initially quasi-two-dimensional electron is becoming increasingly quasi-zero-dimensional. We thus 
expect that the binding energy will increase substantially in this situation. On the other hand, if the magnetic field is 
aligned in the x-direction, magnetic localization involves the y- and z-directions with characteristic lengths A c . Since 
the QW also confines the electron motion in the z-direction, and does not restrict its motion along x, the combined 
effects of magnetic and of QW localization now result in a one-dimensional motion. As seen in Fig. even up to 
73=16 T the magnetic localization now has practically no effect, and binding energy practically does not depend on 
B. In contrast with L=100 A, for a wider QW (L=300 A) the changes in binding energy produced by the magnetic 
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FIG. 3: Donor binding energy E b in a CdTe/Cdo.8iMgo.i9Te QW (L=100 A) as a function of a magnetic fields B for 9=0° 
(Case I), 30°, 45°, 60° and 90° (Case 77). In the Case 77, the binding energy is practically independent on the magnetic field 
B {B <16 T). 

field are quite visible at 0=90°, as seen in Fig. [21 For such a wide QW, the change of E b as B increases from to 
16 T is 0.3 7?y*=4 meV. Now L=300 A and the binding energy is only 1.2 Ry* (without magnetic field) so that the 
system is nearly three-dimensional (as << L), in contrast to the two-dimensional character obtained for L=100 A. 
Thus for all values of 8 the magnetic field effectively localizes the particle in both dimensions perpendicular to the 
direction of the applied field. This explains why the curves for E b are much more flat for L=300 A than for L=100 A, 
particularly at higher values of B. We expect that, when we increase the QW width even more, E b should become 
even more flat, eventually approaching the three-dimensional limit, where it ceases to depend on the tilt angle even 
for large B. Our calculations clearly confirm this trend. 

In Fig. [3 we show the dependance of the donor binding energy E b for a CdTe/Cdo.siMgo.igTe quantum well 
(£=100 A) as a function of the magnetic field 73, for five different values of the tilt angle 9. We have chosen the 
magnetic field range to be < B < 16 T, which is the most widely accessible field range in photoluminescence (PL) 
spectroscopy. As seen in the figure, the binding energy at 73=0 T is E b {B ~ 0T) ~ 1.7 Ry*, which corresponds to a 
70% increase in binding energy with respect to the binding energy of a three-dimensional donor. This increase is due 
to the confinement generated by the quantum well. For Case 7, the increase in magnetic field has a clearly visible 
impact on E b (B). In contrast, E b is practically constant for Case 77. While the series of curves presented in Fig. 
seem to be linear, our results for the binding energy have, in fact, a square-root dependence on the external magnetic 
field. Such a dependence is in accordance with the results obtained in Ref. Q. The apparent linearity of the curves is 
due to the fact that at the highest field we consider (73=16 T) the ratio j2Ry* is only 0.74. The E b (B) cx \f~B scaling 
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behavior becomes apparent for a much wider range of magnetic fields < 7 < 5. In Fig. 0] we present results for 
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FIG. 4: Donor binding energy E), in a rectangular CdTe/Cdo.98Mgo.o2Te quantum well as a function of tilt angle 9 of external 
magnetic field B for different magnetic field values (different symbols), calculated for two QW widths L. Full lines correspond 
to L=100 A, and dashed lines to L=300 A. Lines without symbols correspond to B=0 T. 



14=20 meV (corresponding toisi 0.02). Comparing Fig. [21and Fig. 01 we see that the binding energy of the donor 
is larger for 14=200 meV than for 14=20 meV. This well known fact originates from the larger quantum confinement 
of the deeper QW. The characteristics of the results in Fig. 0] are similar to those in Fig. [3 including the monotonic 
behavior of Ef,. Comparing the curves in Fig. [3 with corresponding curves in Fig.0] we see that the latter clearly are 
more flat. This again confirms that in the three-dimensional case, i.e., as V e — * 0, we should have no dependance 
(straight horizontal lines). In Fig. we show the difference Ediff between the binding energy of a donor at 9 = 0° 
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FIG. 5: Difference Ediff of the binding energy of a donor for 9—0 and 0=90° in a quantum well with L — 50 A as a function 
of barrier height Ve. Full symbols correspond to B—7 T, and open symbols to _B=16 T. 

and its binding energy at 6 = 90° as a function of the height of the barrier V e ; i.e., Ediff — E(6 — 0°) — E(6 = 90°). 
At B=7 T and V e > 50 meV, the difference Ediff is practically constant and relatively small (only 3.7 meV), but 
at B=16 T it does not saturate until V e «75 meV; and its value is twice as high, i.e. 7.5 meV. This feature can be 
utilized as a tool for determining the conduction band offset, at least in QWs with moderate barrier heights [22l |. 
Note that the bigger the magnetic field, the larger the offset which can be measured using this method. 
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In Ref. [16(, Fig. 6, the binding energy of the donor as a function of tilt angle is a non-monotonic function showing 
a maximum at 9 = 45°, in contrast to the monotonic behavior reported here. In our opinion this is related to the 
approach employed by the authors of R ef. Ilq , in which a real quantum well is transformed into two QWs oriented 
at right angles to one another. In Ref. |14| the same group, using the same approximation, calculated the exciton 
binding energy as a function of tilt angle (see Fig. 7 in Ref. 14] ). Unfortunately, the approach used in the latter 
reference does not provide the results for the range of 9 between 0° and 15°, and between 75° and 90°, which appears 
to be an artefact of the technique used in Ref. an d R- e f ■ ■ 

IV. CONCLUSIONS 

We have shown the results of variational calculations of the binding energy of a neutral donor in a rectangular 
quantum well as a function of the angle of en external magnetic field tilted with respect to the growth direction of the 
QW. In our calculations we used parameters characteristic of II- VI compounds (using specifically parameters for the 
CdTe/CdMgTe QW system), and we assumed that the donor is located in the center of the quantum well. For a given 
magnetic field, the largest binding energy is found to correspond to the case when the magnetic field is perpendicular 
to the plane of the QW. We find that the binding energy of D° is a monotonic function of the tilt angle 9, decreasing 
with increasing tilt angle, in contrast with earlier calculations reported in Refs. [[Til ITt|]- Our results reduce to the 
three-dimensional limit when either the quantum well width increases or the barrier height decreases, providing a 
"reality check" for the method used. We have shown that for the CdTe/Cdi-^Mg^Te quantum well (x < 0.1), the 
conduction band offset can be determined by measuring the binding energy of the neutral donor at two perpendicular 
directions of the applied magnetic field, 9=0 and 90°. To our knowledge this technique of determining conduction 
band offsets has not been previously recognized. 
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